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ABSTRACT 

We studied the impact of the revisited values for the LSR circular velocity of the Milky Way 
(Reid et al. 2004) on the formation of the Magellanic Stream. The LSR circular velocity was 
varied within its observational uncertainties as a free parameter of the interaction between the 
Large (LMC) and the Small (SMC) Magellanic Clouds and the Galaxy. We have shown that 
the large-scale morphology and kinematics of the Magellanic Stream may be reproduced as tidal 
features, assuming the recent values of the proper motions of the Magellanic Clouds (Kallivayalil 
et al. 2006). Automated exploration of the entire parameter space for the interaction was 
performed to identify all parameter combinations that allow for modeling the Magellanic Stream. 
Satisfactory models exist for the dynamical mass of the Milky Way within a wide range of 
0.6 • 10 12 Mq to 3.0 • 10 12 Mq and over the entire 1-a errors of the proper motions of the Clouds. 
However, the successful models share a common interaction scenario. The Magellanic Clouds are 
satellites of the Milky Way, and in all cases two close LMC-SMC encounters occurred within 
the last 4Gyr at t < — 2.5Gyr and t ss — 150 Myr, triggering the formation of the Stream and 
of the Magellanic Bridge, respectively. The latter encounter is encoded in the observed proper 
motions and inevitable in any model of the interaction. We conclude that the tidal origin of the 
Magellanic Stream implies the introduced LMC/SMC orbital history, unless the parameters of 
the interaction are revised substantially. 



Subject headings: galaxies: evolution — galaxies: interactions 
— Magellanic Clouds — methods: numerical 



galaxies: kinematics and dynamics 



1. Introduction 

Modeling the evolution of the Magellanic 
Clouds and their interaction with the Milky Way 
is a challenging task. We have to deal with a 
unique system of nearby galaxies which have been 
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subject to various detailed observational surveys. 
The resulting amount of data has established an 
extended set of constraints for every theoretical 
study of the formation and evolution of the Mag- 
ellanic Clouds. On the other hand, namely the 
information about the kinematics and the mor- 
phology of the large-scale structures associated 
with the Clouds have proved useful to reduce sig- 
nificantly the parameter space of the interaction 
between the Large Magellanic Cloud (LMC), the 
Small Magellanic Cloud (SMC) and the Milky 
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Way. 

Undoubtedly, the spatial velocities of the Mag- 
ellanic Clouds are the most critical input parame- 



ters f or models of their interaction (jRuzicka et al 
2009). Regarding the difficulties accompanying 



the observational measurements of the proper mo- 
tions of the LMC and the SMC, it is not surprising 
that the results have al ways suffered from large 
uncertainties (see, e .g. iKroupa fe Bastianl Il997 ; 



Pedreros et al.l 120021) . Thus, the most important 



parameters were also the most uncertain ones. 

A common feature of models of the Magellanic 
Clouds (see Section^) are bound orbits around the 
Galaxy, with several revolutions over the Hubble 
time. H owever, recent proper mo tion measure- 
ments by iKallivavalil et al. ( 2006bl la ) have intro- 
duced entirely surprising results. Besides a sub- 
stantial reduction of the errors, they claim sig- 
nificantly increased galactocen tric velocitie s . The 
subsequent detailed analysis bv lBesla et al. ( 2007 ) 
has revealed a new value of the LMC galactocen- 
tric velocity Wlmc ~ 350 km s -1 , which is close to 
the local escape velocity. Thus, the scenario with 
several perigalactic approaches has become rather 
unlikely. 

Through a detaile d automated explorat ion of 
the parameter space, iRuzicka et ail (l2009h have 



shown that the new data by IKallivavalil et al 



(|2006bl lah rule out tidal stripping as the process 
possibly dominating the formation of the Magel- 
lanic Stream - a gaseous trailing tail emanating 
from the Magellanic Clouds, crossing the South 
Galactic Pole and stretched oyer ~ f00° of the 
plane of sky. iMastropietr 3 (|2009l ) later claimed 
that the efficiency of ram pressure increased due 
to the higher relative velocity of the Clouds and 
the hot halo gas. This might compensate for the 
shorter timescale for the interaction, and allow for 
the creation of the Magellanic Stream of a spatial 
extent comparable to the observations. 

The analysis of the orb ital history of the LMC 
bv lShattow fc Loebl (|2009h allowing for the modi- 
fied values of the LSR circular velocity 6n and the 
LSR angular rotation rate Slo (|Reid fc Brunthaler 
2004) has introduced promising results. The com- 
bination of changes to both, the kinematics and 
dynamics of the interaction, may result in the de- 
crease of the total energy of the LMC-Milky Way 
pair, and keep the LMC within the tidal radius of 
the Galaxy for several Gyrs. 



Following the findings by IShattow fc Loeb 
(120091 ). we have investigated the dynamical evo- 
lution of the Magellanic Clouds and their inter- 
action with the Galaxy in terms of a fully au- 
tomated exploration of the parameter space for 
the interaction. The goal of this study is to show 
whether the up-to-date observational estimates 
of the LMC/SMC proper motions together with 
the revisited values for the solar galactocentric 
motion and distance allow for the reproduction of 
the large-scale distribution of neutral hydrogen 
(HI) in the LMC-SMC-Milky Way system (the 
Magellanic System) due to the gravitational tidal 
strip ping. The approach app lied was introduced 
bv lRuzicka et all (20071 120091) . and it is based on 
the evolutionary optimization (genetic algorithm) 
of the model input according to its ability to re- 
produce the_JJI_oteervations of the Magellanic 
System (jBriins et al.ll2005l ). 

The procedure devised bv lRuzicka et al. I (120091) 
was followed with a couple of modifications. The 
ranges of the proper motions of the Clouds were 
redefine d in order to accommoda te not only the 
data bv IKallivavalil et al.l (2006b. a) but also the 
results by Piatek et al.l ( 20081) who reprocessed in- 
dependently the original raw proper motion data 
(Figure [TJ. The parameter space of the LMC- 
SMC-Milky Way interaction was also extended by 
an additional dimension representing the LSR cir- 
cular velocity. A more detailed description of the 
studied parameter space will be provided in Sec- 
tion [3l 



The read er familiar with the papers by Ru zicka et al 
( 2007 . 20091 ) may skip Section |1] which introduces 
the technical details of our numerical model and 
the principles of the automated exploration of pa- 
rameter spaces. In Section [5] the influence of the 
LSR circular velocity on the interaction between 
the Magellanic Clouds and the Galaxy is stud- 
ied. Section [6] offers our results which arc later 
summarized in Section [7J followed by a discus- 
sion in Section |8l Finally, the most challenging 
open questions related to the formation of the 
large-scale Magellanic structures are mentioned 
in Section [9] 

2. Models of the Magellanic System 

It has been a common practice in the theoreti- 
cal studies of the Magellanic System to adopt sev- 
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eral reasonable assumptions on the orbital paths 
of the Clouds to reduce the volume of the cur- 
rent LMC/SMC velocity space. In most cases the 
Magellanic Stream has been used as a tracer of 
the past orbits o f the M agellanic Cloud s (e.g . 
Fuiimoto fc Sofuel Jl976t iGardiner et~aH Il994t 
Connors et al.ll200i 



This approach has allowed for remarkably good 
models of the large- scale distributi o n of H I in the 
Magellanic System. IConnors et al. I (120061) devised 
a complex tidal model of the interaction between 
the Clouds and the Milky Way offering an impres- 
sive reproduction of the kinematics of HI in the 
Magellanic Stream. An over-density was present 
at the tip of the Stream, however, demonstrating 
the common problem of the tidal stripping. 

Alternatively, ansatz explaining the origin of 
the Stream by the ram pressure stripping of the 
LMC/SMC neutral hydrogen has been introduced 
as an explanation. The adv ance d numerical mod- 
els bvlBekki fc Chibal(|2005l) and lMastropietro et a . 
(2005) simulated the distribution of HI at the 



main body of the LMC and the global morphol 



ogy o f the Magellanic Stream. iBekki fc Chiba 
(2005) described the LMC as a system of self- 
gravitating particles. To account for the dissipa- 
tive processes in the gaseous medium, the method 
of sticky particles was emplo yed. The SPH code 
bv lMastropietro et al.l (|2005h allowed the redistri- 
bution of H I in the Clouds to occur at timescales 
of the order 10 s yr only. 

3. Parameter space of the Magellanic Sys- 
tem 

Various aspects of modeling observed interact- 
ing galaxi es were discus sed to a great deta il in th e 
studies bv lTheid (Il999h or iTheis fc Kohld (l200lh . 
Obviously, even an interaction of two galaxies de- 
scribed in a strongly simplified manner involves 
well over 10 parameters necessary to describe the 
structure, the total mass and the mass distribution 
of both galaxies, and also their initial positions 
and velocities. In principle, the values of the pa- 
rameters may be derived from observational data. 
However, regarding distances of galaxies, acquir- 
ing observational data of the resolution and variety 
sufficient to uniquely determine the parameters of 
their interactions is very difficult. Unfortunately, 
this is true even for the satellites of the Milky Way. 



Thus, modeling interacting galaxies involves in- 
evitably an exploration of the parameter space for 
the interaction to determine the parameter combi- 
nations, i.e. the evolutionary scenarios, compati- 
ble with the observed kinematics and morphology 
of the galactic system. 

3.1. Key parameters 

The number of parameters required to describe 
the interaction involving the Milky Way, the LMC 
and the SMC depends on the physical model 
adopted. We have employed a restricted N-body 
numerical code assuming the gravitational inter- 
action only. Even such a simplified view of the 
interaction involves ~ 25 independent parame- 
ters including the initial conditions of the LMC 
and the SMC motion, their total masses, param- 
eters of mass distributions, particle disk radii, 
and orientation angles, and also the parameters 
defining the gravitational potential of the Milky 
Way, i.e. the flattening parameter of the dark 
matter (DM) halo, the LSR circular velocity and 
the LSR angular rotation rate (see Section [5]). 
Some of the parameters were constrained by the- 
oretical studies (including scale radii e of the 
LMC/SMC halos, the Coulomb logarithm A for 
the dynamical friction in the Milky Way halo, and 
the halo flattening parameter q for the model of 
the gravitational potential of the Galaxy). Their 
m ean values and searc hed errors were discussed 
in IRuzicka et al.l ( 2007 ). This section introduces 



the observationally estimated parameters of the 
LMC-SMC-Galaxy interaction (see also Table [[]). 

It has been shown by IRuzicka et al. (|2009l ) 
that the current proper motions of the Magellanic 
Clouds are the most critical parameters of their in- 
teraction with the Milky Way. Figure [T] shows the 
proper motions of the LMC and the SMC as esti- 



mated by the l atest analysis bv |Kallivavalil et al. 
(|2006bl lari and IPiatek et all (j2008|) based on the 



measurements by the HST. For convenience, the 
LMC/SMC proper motion vectors are decomposed 
into two components related to the local direc- 
tions of west (^w) and north (^n)- The rea- 
son for involving two results derived from a sin- 
gle da ta-set is the discr e pancy betw een the stud- 



i es by iKallivavalil et al.l (|2006a) and IPiatek et al 



(|2008l ) regarding the SMC western proper motion 
/^w- The resulting proper motion space for the 
Magellanic Clouds explored in this study is a union 
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of the sets defined b y the 1- a error s of the data 
bv iKallivavalil et al.l (|2006bl lal) and IPiatek et all 
j2008|) for the LMC and the SMC, respectively 
(Figured]): 



line 
lmc 

smc 
smc 

Mn 



(-2.11, -1.92) mas yr" 1 
(+0.39, +0.49) mas yr" 1 
(-1.34,-0.69) mas yr" 1 
(-1.35,-0.99) mas yr" 1 . 



(1) 



Detailed analysis by IRuzicka et al. (l2009h has 
quantified the sensitivity of the dynamical evolu- 
tion of the Magellanic System to the variation in 
different free parameters. Except the LMC/SMC 
current spatial velocities and positions, the re- 
maining parameters had a rather weak impact 
on the interaction. As long as the parameters 
are independent adding another parameters to 
the original set does no t change the conclusion 
by IRuzicka et al. ( 20091) regarding their impor- 
tance. However, once parameter independence 
cannot be guaranteed, it is a good practice to keep 
the dimensionality of the parameter space as high 
as possible. Therefore, we will introduce the LSR 
circular velocity as a new dimension of the Magel- 
lanic parameter spac e established in the way sim- 
ilar to IRuzicka etail (|2007i [20091) . 

Unlike the proper motion of the Clouds their 
LSR radial velocities could be measured wi t h high 



accuracy. Following Ivan der Marel et al.l ((2002) 



we set v 1 ^ = 262.2 ±3.4 km s" 1 . The SMC radial 
veloci ty error was estimated by lHarris fe Zaritskv 



(2006) 



as v 



rail 



146.0 + 0.60 km s" 



was adopted from Ivan der Marel et al 



The heliocentric position vector of the LMC 

(l2002h . 

i.e. the equatorial coordinates are (oii mc , <5i mc ) = 
(81.90°±0.98°, -69.87°±0.41°), its distance mod- 
ulus is (m-M) lmc = 18.5 + 0.1. The equa- 
torial coordinates of the SMC were set to the 



ranges (a E 



(13.2° ± 0. 3°, -72.5° ± 
0.3°) (see IStanimirovic et al.l |2004 and refer- 



si ic i "smc I 



ences therein). A range of dista nce determina- 
tions f or the SMC was provided bv lvan den Bergh 
(2000) and we used his resulting distance modulus 
(™-M) smc = 18.85 ±0.10. 

Several observational determinations of the 
LMC disk pla ne orientation have been published 
so far (see, e.g. Lin et al . 1995: Ivan der Marel et al. 
20021) . In our parameter study the LMC in- 
clination i and position angle p and their er- 



rors agree with Ivan der Marel et al.l (|2002j), i.e. 
i = 34.7° ± 6.2° and p = 129.9° ± 6.0°. As 
the SMC misses a well defined disk, the orien- 
tation and the position an gle usually refer to 
t he S MC "bar" defined by iGardiner fc Noguchi 
1996). B ased on the estimat e s by yan den Bergh 



20001 ) or IStanimirovic et al.l (|2004l) we adopted 
the error ranges i = 60° ± 20° and p = 45° ± 20° 
for the SMC initial disk inclination and position 
angle, respectively. 

Gardiner et al.l ( 19941) analyzed the HI surface 



contour map of the Clouds to estimate the initial 



LMC and SMC disk radii 



disk 



and 



r disk> respec- 



tiv ely. With the use of their work and of the study 
bv iBruns et ail (|2005l) we varied the LMC/SMC 
disk radii within the ranges of 10.5±1.5 kpc (LMC) 
and 6.5 ± 1.5 kpc (see Table [1]). Regarding the ab- 
sence of a clearly defined disk of the SMC and pos- 
sible significant mass redistribution in the Clouds 
during their evolution, the results require careful 
treatment. 

Current t otal masses m\ mc and m smc follow the 
estimates bv lvan den Berghl(l2000l) . The masses of 



the Clouds are functions of time and evolve due to 
the LMC-SMC exchange of matter and as a conse- 
quence of the interaction between the Clouds and 
the MW. Our test-particle model does not allow 
for a reasonable treatment of a time-dependent 
mass loss. Therefore, the masses of the Clouds 
are considered constant in time and their initial 
values at the starting epoch of simulations are ap- 
proximated by the current LMC and SMC masses. 

The gravitational potential of the Milky Way 
is modeled by the superposition of three static 
components includi ng the axially symmetric l og- 
arithmic potential (|Binnev fe Tremaind 1987 ) of 
the DM halo, the Miyamoto-Nagai potential of 
the Galactic disk, and the Hernquist bulge. Since 
the logarithmic halo is the dominant component 
of our model, its flattening parameter was treated 
as a free parameter. The flattening q was var- 
ied within the range of (0.71,1.30). The lower 
limit is based on the condition q > 1/ \/2 required 
to avoid negative mass densities in the logarith- 
mic halo. The upper limit was set up to intro- 
duce a convenient symmetry with respect to the 
spherical shape (q = 1). The scaling velocity fac- 
tor vq of the logarithmic model (more details in 
Binnev fc Tremainelll987 ) is a function of the LSR 



circular velocity and the solar galactocentric dis- 
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tance. 

The freedom in the values of the parameters q 
and vo introduced a spread in the total mass of the 
Galaxy m M w = (0.59,5.90) • 10 12 M Q within the 
radius of 250 kpc. The estimates of the total mass 
of the Milky Way exceeding sa 3 • 10 12 M Q have 
neit her observational not theoretical support (se e 
e.g. iBinnev fc Tremaind [l987t iLi fe Whitd 120081) . 
but this fact will be taken into account when the 
results are interpreted. 

3.2. Rotation curve of the Milky Way 



Unlike the studies by iRuzicka et al. ( 2007 



2009t ) we have treated the LSR circular velocity 
0o as a free parameter. Thi s decision was m o- 



tivate d by recent findings by iReid fe Brunthaler 
([20041 ) . Their measurements of the proper motion 
of Sag A* yielded the estimate for the LSR circular 
velocity 8 = 236 ± lSkms -1 . At the same time, 
they give a revisited value of the LSR angular ro- 
tation rate fio of 29.45 kms -1 kpc -1 . These two 
quantities then define the variations of the solar 
galactocentric distance, as Rq = Qq/^o- 

IReid fe Brunthalerl (|2004) performed very pre- 
cise position measurements of the famous radio 
source Sag A* with respect to two background ex- 
tragalactic radio sources. The data collected over 
a period of 8 years yielded values of the appar- 
ent motion of Sag A*. This quantity is a com- 
position of the reflected solar galactocentric mo- 
tion and of the peculiar motion of Sag A* which 
is quite small. The fraction of the apparent mo- 
tion corresponding to Solar rotation around the 
Galactic Center may be treated as a combination 
of the Local Standard of Rest (LSR) circular ve- 
locity Go and the devi ation of the solar motio n 
from the circular orbit dDehnen fe Binnevl ll998h. 
Finally. [Reid fe Brunthalerl (|2004f ) were able to de- 
termine the latter which yielded an estimate for 
the LSR circular velocity Oo = 236 ± 15kms _1 
assuming the solar distance to the Galactic cen- 
ter (i?o) of 8.0 ±0.5 kpc. Notably, the uncertainty 
of Qq originating from the measurement error was 
only lkms -1 . The rest appears due to the un- 
certainty of i?o (|Reid fe Brunthalerll2004l) . Calcu- 
lating the LSR angular rotation rate f^o = ®o / Ro 
yields a value with a very small uncertainty: Qq = 
29.45 ± 0.15kms _1 kpc -1 . The corresponding rel- 
ative error is only 0.15/29.45 = 5 • 10~ 3 which 
is by one order of magnitude below the relative 



error of any of the parameters examined in our 
study. Therefore, the LSR angular rotation rate 
fio will be treated as a constant in this paper, i.e. 
tt = 29.45 kms" 1 kpc" 1 . 

Although we are already familiar with the 
impact of most of the discussed parameters on 
the behavior of the restricted N-body model 
of the LMC-SMC- Milky Way interaction (see 
Ruzicka et alj|2009h . this does not apply on the 
LSR circular velocity. Even though it is a nat- 
ura l choice, taking t h e 1-a error box of 0o 
by IReid fe Brunthalerl (|2004f ) for the range ex- 
plored by the parameter study is arguably not a 
good practice. We want to avoid the case of im- 
portant models being localized at either end of 
the studied parameter range. Such a coincidence 
would then require an extension of the studied 
range of the LSR circular velocity followed by 
another time-consuming run of the automated 
parameter search. The described procedure is 
necessary to uncover and understand the behavior 
of the interacting system in such a parameter re- 
gion of apparently special features. Therefore, the 
range of the new input parameter Oo was obtained 
by extending the 1-a uncertai nty of the LSR cir- 
cular velocity as published in Reid fc; Brunthalerl 



(|2004fl by lOkms" 1 , i.e. 9 = (210, 260) km s" 



corresponding to the range of the Solar galacto- 
centric distance i?o of (7.13, 8.83) kpc. 

4. Method 

This paper focuses on the gravitational interac- 
tion between the Magellanic Clouds and the Milky 
Way. Although hardly any doubts exist regard- 
ing the presence of the hydrodynamical processes, 
as the LMC/SMC neutral hydrogen interacts with 
the ambient hot gas of the Galactic halo, the tidal 
stripping is intrinsically involved in every model 
assuming the gravitational interaction. 

4.1. Restricted N body model 

The model itself is an advanced version of 
the scheme by iGardiner et al.l ( 19941 ) : it is a re- 
stricted N-body (i.e. test particle) code describing 
the gravitational interaction between the Galaxy 
and its dwarf companions. The potential of the 
Milky Way is dominated by the flattened dark 
matter halo and dynamical friction is exerted on 
the Magellanic Clouds as they move through the 
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halo (|Binnevl [l977h . The LMC and the SMC 

are represented by Plummer spheres, initially sur- 
rou nded by test-particle disks . For further details 



see 



Ruzicka et atl (|2007l [2009). 



4.2. Searching the parameter space 

The exploration of the parameter space for the 
interaction involving the Galaxy and the Mag- 
ellanic Clouds was performed by a genetic algo- 
rithm. Genetic algorithms belong to the class 
of evolutionary optimizers that mimic the selec- 
tion strategy of natural evolution. iHollandl (1975) 
was the first one who proposed the application 
of genetic algorithms on optimization problems 
in mathematics. Recently, the performance of 
genetic algorithms was studied for galaxies in in- 
teract ion (see, e.g. iTheid fl999l) . iTheis fc KohTel 
(|200lh analyzed the parameter space of two 
observed interacting galaxies - NGC4449 and 
DDO 125. Genetic algorithms turned out to be 
very robust tools for such a task if the rou- 
tine comparing the observational and modeled 
da ta is appropria t ely d efined. The approach 
bv lTheis fc Kohld (l200lh was later adopted and 
improved in order to explore the int eraction of the 



Magellanic Clouds and the Galaxy (IRuzicka et al 
2007ll2009[ ). 



The comparison between the model and obser- 
vations became more efficient by invo lving an ex- 
plicit search for the shapes in the data ([Ruzicka et al 
2009J). Also the significant system-specific fea- 
tures (such as a special geometry and kinemat- 
ics) were taken into account, further improving 
the performance of the genetic algorithm for ex- 
ploration of the LMC-SMC-Galaxy interaction. 
More detailed information is to be found in Sec- 
tion H3] of this paper. 



4.3. Fitness function 

The automated search of the parameter space 
is driven by a routine comparing the modeled 
and observed distributio ns of H I associate d with 
the Magellanic System (jBriins et al.ll2005l) . The 
match is measured by the fitness function (/) 
which is a function of all input parameters, as ev- 
ery parameter set determines the resulting simu- 
lated H I data-cube. Our fitness function returns a 
floating-point number between 0.0 (complete dis- 
agreement) and 1.0 (perfect match) and consists 



of four different comparisons, including search for 
structures and analysis of local kinematics. For 
details on the fitness function used for this study 
see Appendix |A1 

The fitness function / is the only part of the 
genetic algorithm that reflects the nature of the 
studied problem. Although the automated evolu- 
tionary optimization is a robust method applicable 
on a remarkable variety of systems, its actual per- 
formance and efficiency are critically dependent 
on the proper choice for the fitness function. If 
the function is defined in a sensible way, the num- 
ber of parameter combinations examined by the 
genetic algorithm is minimized. However, the per- 
formance of the genetic algorithm also depends on 
its convergence rate expressed as df(i)/di, where 
i = 1,2,... stands for the number of the actual 
generation and f(i) is the mean value of the fit- 
ness function in the i-th generation. The conver- 
gence rate is generally proportional to the ratio of 
the generation size JV gcn and of the dimension n 
of the studied parameter space over which the fit- 
ness function is defined. Unfortunately, the value 
of N gon /n is very low if genetic algorithm-based 
optimizer is applied on the parameter spaces of in- 
teracting galaxies. The number of parameters in- 
volved is always high (see Section[3]). On the other 
hand, the maximum generation size has to be lim- 
ited in order to keep the total computational time 
requested by the parameter search reasonable. 

In order to overcome such a difficulty we 
searched the parameter space repeatedly in a fixed 
number of optimization steps, i.e. generations of 
models, and localized « 100 high-fidelity models. 
Such a procedure is not likely to reveal the global 
maximum of the fitness function, but it results in a 
map of distribution of good models over the entire 
parameter space. In principle, every region of the 
parameter space allowing for the satisfactory re- 
production of the large-scale structures associated 
with the Magellanic Clouds may be identified. 

The reader might ask why there is such em- 
phasis placed on the fitness function itself if it, in 
fact, does not seem to provide any physical infor- 
mation about the interacting system of the Galaxy 
and the Magellanic Clouds. Indeed, the function 
/ serves primarily as a driver to the genetic algo- 
rithm. However, the search for good models of the 
observed Magellanic System is efficient only if rel- 
evant astrophysical data are supplied as the input 



G 



to /. As already mentioned, our study deals with 
detailed morphological an d kinematic inform ation 
from the 21 cm survey bv lBriins et all (j2005h and 
with the corresponding modeled data. The fitness 
function then makes a link between the observ- 
able data and the initial state of the Magellanic 
System. 

5. The role of the Solar galactocentric ve- 
locity and distance 



The conclusions by iReid fc Brunthalerl (j2004l) 
regarding the LSR angular rotation rate f2 and 
the LSR circular velocity O have a significant im- 
pact on both the kinematics and the dynamics of 
the Magellanic System. 

The dynamics of the interaction between the 
Magellanic Clouds and the Galaxy is influenced 
by the change of the mass distribution and of the 
to tal mass of the Milky W ay. The value of CIq 
bv lReid fc Brunthalerl (|2004l) implies that the IAU 
standards for 0q AU and for the Solar galactocen- 
tric distance i?g AU cannot hold at the same time 
anymore, as 

220kms~ 1 , , 

= 25.88 km s^kpc" 1 



o IAU 



Rl AV 



8.5 kpc 



while IReid fc Brunthalerl (|2004l ) expect that flo = 
29.45 kms -1 kpc -1 . Therefore, a rescaling of the 
rot ation curve of the Galaxy occurs if the results 
by IReid fe Brunthalerl (2004) are taken into ac- 
count, which necessarily means a rescaling of the 
mass distribution in the Galaxy. 

The current positions of the Magellanic Clouds 
in the phase space have been measured with re- 
spect to the phase space position of the Sun. In 
order to solve the equations of motion for the inter- 
action between the Clouds and the Milky Way, the 
galactocentric positions and velocities of the LMC 
and the SMC are needed, i.e. the knowledge in the 
LSR circular velocity 0o and the LSR galactocen- 
tric distance Rq is necessary (s ee Figure Et . Thus, 
i t is ob vious that the results bv lReid fc Brunthaler 
([20041 ) may influence the kinematics of the Mag- 
ellanic Clouds. 



5.1. Model of the Milky Way 

The rotation curve of a galaxy 
d$ tot (R,z) 



v 2 (R) = R 



OR 



(2) 



introduces a fundamental constraint on the models 
of its overall potential $tot(-R, z). In our study, the 
following three component model of the potential 
of the Milky Way has been used: 

$ tot (i?, z) = $ L (i?, z) + $ MN (i?, z) + $bM, (3) 

where 

$ h (R,z) = ±vlln(Rl+R 2 + ^\ (4) 
is the logarithmic model of the Galactic DM halo, 



GMa 



\J r 2 + (b + Vz^Tc^y 



(5) 



is the Miyamoto-Nagai potential of the Galactic 
disk with the total mass Md, and 



$bW = -G 



(6) 



is the Hernquist formula for the spherically sym- 
metric potential of the Milky Way bulge with the 
mass Mb- 

At the galactocentric distance where the Mag- 
ellanic Clouds have resided, i.e. D <; 50 kpc, the 
gravitational field of the Galaxy is dominated by 
its DM halo component. The role of the DM halo 
becomes more prominent with the increasing dis- 
tance to the galactic center, because V<I>l ~ 1/r 
while V$mn ~ and also V$b ~ l/r 2 . The 
necessary adjustment of the Milky Way potential 
^tot required by the Equation ([2]), if the LSR cir- 
cular velocity 0o is changed, was achieved by vary- 
ing the velocity constant u in Equation (f?)). 

The choice for the remaining p arameters in the 



Equations dU), ©, and Q follows iFellhauer et al 



z=0 



(|2006l ) and iRuzicka et alj (|2007l 120091 ) and so we 
have set a = 0.7 kpc, b = 6.5 kpc, c = 0.26 kpc, 
M b = 3.41O 1O M , M d = 10 11 M©, and R c = 
12.0kpc. 

The value of the characteristic radius R c de- 
termining the concentration of the Galactic halo 
is five times smaller than the present galactocen- 
tric distances of the Magellanic Clouds, and by 
the order of 10 smaller compared to the typical 
LMC/SMC galactocentric distances over the last 
several Gyr. Hence, its impact on the density pro- 
file of the region of the DM halo of the Milky Way, 
where the Clouds have resided, is weak. Regard- 
ing the total mass of the Galaxy enclosed within 
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a given radius, it is determined by the velocity 
parameter vq, as 



/ PL (R,z,q)d 3 r + M b 

J halo 



(7) 



where the density of the logarithmic halo of the 
Galaxy ph(R, z) ~ v%, while ph(R, z) ~ l/i?c- We 
did not consider the scale radius i? c as a free pa- 
rameter because its influence on the mass distri- 
bution and the total mass of the DM halo of the 
Milky Way is significantly lower compared to the 
velocity parameter vq. 

The flattening q of the DM halo potential en- 
tered our simulations of the Magellanic System as 
a free parameter (see Sec. [3]). The velocity param- 
eter vq of the logarithmic halo was considered a 
function of the LSR circular velocity Qq. The cor- 
responding formula comes out of the Equations ([2]) 
and ([3]) after substituting for $l, ^mn, and $b: 



e i 



Rl 



G [R, + i?n 



M h 



(R 2 + (b+cry 



R (R + a)2 



(8) 



where Rq = 0o/f2o. 

The analysis of the Equation © shows that if 



a < 6q/^o 



(9) 



and 



{b+c-VU>R c ){b+c+VTl>R c ) < 0.59o/^o ( 10 ) 



then 



9vp 



> o. 



It can be easily seen that the conditions (|9]) 
and (Tit)]) are always satisfied for the values of 
the parameters a, b, c, and R c assumed in our 
study. As we have Equation for the loga- 
rithmic halo ( Binnev &: Tremainelll987 ). the total 
mass of the Milky Way enclosed within an ar- 
bitrary radius r must grow if the LSR circular 
velocity 0o is increased. 

5.2. Galactocentric positions and veloci- 
ties of the Magellanic Clouds 

A convenient transformation of the position and 
velocity vectors of the Magellanic Clouds from the 



heliocentric spherical coordinates to the galacto- 
cen tric Cartesian frame (Figu re [2J was derived 



by Ivan der Marel et all (|2002h . Let (l,b) and 
(a, 5) be the corresponding pair of the galactic 
and equatorial coordinates of the center of mass 
of either of the Magellanic Clouds, respectively. If 
the Cloud's present heliocentric distance is D, its 
galactocentric Cartesian coordinates are 



0,1,2, 



(11) 



where Rq is the Solar galactocentric position vec- 
tor and Uo is the heliocentric unit position vector 
of the Cloud. 

The heliocentric spatial velocity of an object in 
space is usually expressed in terms of its proper 
motion and the line-of-sight systemic velocity 
The proper motions in the directions of 



Jsys- 



west and north are defined as ( van der Marel et al 



2002) 



dd 
dt 



37- (12) 



The transformation from the heliocentric velocity 
coordinates /xw, A*Nj and v sys to the galactocentric 
velocity components may be then expressed as 

v l = Vq + v sys UQ + Dp w u\ + Dp, N ul, (13) 

where 



u = 

Ul = 

u 2 = 



(cos I cos b, sin I cos b, sin b) 
1 du 
cos (5 da 

dug 

85 ' 



(14) 



and vq is the Solar galactocentric velocity vector. 

Taking Equations ([TTj) and ([13")) into account 
together with the fact that R© = (— 8o/^o,0, 0) 
and vq — (O,0o,O) (see Figure [2]), it is obvious 
that the current galactocentric positions and ve- 
locities of the Magellanic Clouds depend on the 
values of the LSR circular velocity Qq and of the 
LSR angular rotation rate £Iq. If either of these 
two parameters is varied, a change to the galacto- 
centric position vectors of the Clouds occurs due 
to the change of their ^-components: 



-6oM 



Du° . 



(15) 



Similarly, the galactocentric spatial velocities of 
the LMC and the SMC would be influenced, be- 
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cause their y-components involve the explicit de- 
pendence on Oo: 

v 1 — O + v sys ul + Dfi w u{ + D^u\. (16) 

Note that the phase space configuration of the 
system involving the Sun, the Magellanic Clouds 
and the Galactic Center is such that v 1 ~ 0o, 
while Iv 1 ] ~ — 0o- As a consequence, the magni- 
tude of the galactocentric velocity 



(6 + v sys ul + Dfi w u\ + D^u\) 2 (17) 



of the LMC decreases if the LSR circular velocity 
is increased. Figure [3] reveals clearly that the same 
conclusion holds for the SMC. 

Variation of the LSR circular velocity O and 
of the LSR angular rotation rate S7q modifies both 
the magnitudes and the directions of the position 
and velocity vectors of the Clouds. Specifically, 
increase in the LSR circular velocity causes the re- 
duction of the magnitude of the LMC/SMC galac- 
tocentric velocities. 

6. Results 



The complex nature of the interaction involving 
the Magellanic Clouds and the Milky Way results 
in a high-dimensional parameter space. We have 
studied the parameter space by employing genetic 
algorithms as robust optimization tools character- 
ized by reliability and low se nsitivity to local ex- 
tremes ( Theis fc Kohlei 2001 ). This method was 
adopted to analyze the performance of pure tidal 
models for the Magellanic System. 

Regarding the insufficient convergence rate of 
our genetic algorithm (see Section I4.3|) , the au- 
tomated exploration of the parameter space was 
performed repeatedly and the properties of the 
resulting high-fidelity models were analyzed sta- 
tistically. We have identified the regions of the 
parameter space where promising models of the 
LMC-SMC-Milky Way interaction exist. 

This paper focuses on the influence of the LSR 
circular velocity Qq on the formation of the Mag- 
ellanic Stream provided the current spatial veloc- 
ities of the Magellani c Clouds respect the rece nt 
proper motion data bvlKallivavalil et al. I (l2006bllah 
and bv lPiatek et al. ( 20081) . In particular, we want 
to address the question of how to choose this pa- 
rameter in a way that agrees with the recent mea- 



surements by iReid k, Brunthaler ( 2004 ) and al- 
lows for the reproduction of the observed kine- 
matics and spatial extent of H I in the Magellanic 
Stream. 

6.1. Exploration of the parameter space 

With the use of the genetic algorithm, ~ 10 6 
parameter combinations, i.e. individual restricted 
N-body simulations, were examined in total, and 
~ 100 sets providing the highest fitness models 
were collected. In this section the features of these 
models will be discussed with respect to the pa- 
rameters influencing the motion of the Magellanic 
Clouds. 

We have treated the LSR circular velocity as 
a free parameter in order to analyze its impact 
on the tidal interaction in the Magellanic Sys- 
tem. The automated exploration of the modi- 
fied parameter space revealed a significant qual- 
itative change regarding the preferred proper 
m otions of the Clouds , compared to the work 
bv IRuzicka et all (|2009h . 

They have done a detailed analysis of how 
the models of the LMC-SMC-Galaxy interaction 
depend on different parameters. It was clearly 
shown that the interacting system is very sen- 
sitive to the parameters defining the past posi- 
tions of the Clouds' orbits in ph ase space. In 
the work by IRuzicka et al. ( 20091) the set of the 
critical parameters involved the current heliocen- 
tric positions and velocities of the Magellanic 
Clouds and, to some extent, the flattening pa- 
rameter of the Milky Way halo. However, intro- 
ducing t he new value of the L SR an gular rotation 
rate by Reid fc Brunthaler] ( 2004 ) and treating 
the LSR circular velocity as a free parameter has 
changed the old picture. 

Figure |4] shows the distribution of the local 
peaks of the fitness function / (i.e. the fitness of 
every model identified by the genetic algorithm) 
over the studied ranges of the western and north- 
ern components of the proper mo tions for both 
Magel lanic Clouds. In contrast to IRuzicka et al 



(|2009h . ;here is no apparent preference for spe- 
cific values of the proper motions for either of the 
Clouds. Models of a very similar quality have 
been located over the entire studied proper mo- 
tion ranges of the LMC and the SMC. Notably, the 
high-fitness models fall into either of two groups 
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clearly separated by the fitness values of their 
members. The majority of the models has the fit- 
ness exceeding the value of 0.55 while most of the 
rest are described by the relation / < 0.50. Only 
~ 5% of all the models are located within the fit- 
ness range of 0.50 < / < 0.55. Later in Section [B~7l 
the lower limit for the fitness of the satisfactory 
models will be precisely established and discussed. 
For the moment we will just mention that the al- 
most empty belt between the fitness values of 0.50 
and 0.55 separates the acceptable models of the 
Magellanic System fro m those reproducing the H I 



(2005) insuffi- 



observational data by iBriins et al 
ciently. 



6.2. LSR circular velocity and the halo 
flattening as parameters 

It can be easily seen from Equations dU and flSJ) 
that the LSR circular velocity Go and the flatten- 
ing q of the Galactic DM halo determine the to- 
tal mass of the Milky Way. We have focused on 
the distribution of all the high-fidelity models of 
the LMC-SMC-Galaxy interaction over the q— Go 
plane of the studied parameter space. Figure [S] vi- 
sualizes the result. The actual value of the total 
enclosed mass of the Milky Way does not intro- 
duce any limitation regarding the ability of our 
tidal model to reproduce the observed HI large- 
scale Magellanic structures. Such behavior exists 
due to the fact that local peaks of the fitness func- 
tion (i.e. good models) of a similar quality have 
been localized over the entire range of the LSR 
circular velocity. 

The distribution of the genetic algorithm fits 
with respect to the flattening parameter q is quite 
different. A strong preference for oblate (q < 1) 
configurations of the gravitational potential of the 
Milky Way halo ex ists. Such a re s ult ag rees with 
the conclusions bv iRuzicka et al. ( 2007 ). In the 
next paragraphs the role of the LSR circular ve- 
locity and the flattening of the Galactic halo will 
be addressed regarding the impact of these pa- 
rameters on the orbital history of the Magellanic 
Clouds. 



6.3. Galactocentric velocity of the LMC 
and the mass of the Milky Way 

We have revealed an interesting anti-correlation 
in the output of the automated parameter search 



by the genetic algorithm. Figure [6] depicts the 
relation between the magnitude of the current 
galactocentric velocity of the LMC and the total 
mass of the Milky Way for all the high-fitness 
models. The general trend of the dependence, i.e. 
the galactocentric velocity of the LMC decreasing 
as the Galactic mass increases, is driven by the 
LSR circular velocity Go due to the way it links 
the Equations (O, (j8]) and (fl"7|) . and together with 
the fact that the mass of the Milky Way is domi- 
nated by the DM halo. 

Equation ([7|) may be expressed as 

muw = v$(Q )g(q) + M h + M d , (18) 

and for the galactocentric velocity we get 



h(@ , a- u (Sj, (to - M) u fAff, /4, v l iad ), (19) 



where i = lmc. Equations (1181) and (|T9|) define 
the parametric representation of a curve in the 
wmw — |uimc | plane. If all parameters except the 
LSR circular velocity Go are assumed constant and 
Go is varied within the range of (210, 260) kms -1 , 
the parametric equations yield the curves depicted 
in Figure |6] The curves are associated with ran- 
domly selected high-fitness models by taking their 
parameter values for the input to Equations ()18j) 
and ([TO)) . Apparently, the curves are close to lin- 
ear. Their slope is controlled by the LSR circular 
velocity and they are well aligned with the dis- 
tribution of models in the tomw — |^imc| plane 
which provides the key to the Figure H3 The dis- 
tribution of models in the plot is primarily con- 
trolled by the LSR circular velocity which sets the 
linearly decreasing trend of the velocity-mass de- 
pendence. The vertical (velocity) and horizontal 
(mass) placement of the curves depends on the he- 
liocentric positions and velocities of the LMC and 
on the flattening of the DM halo of the Milky Way 
respectively. 

6.4. Orbits of the Magellanic Clouds 

Figure [7] will help us to understand how the 
variations of the flattening parameter q and of the 
LSR circular velocity are reflected in the past or- 
bits of the Clouds. One of the parameters was 
varied within its entire range considered in this 
study, while the second one was fixed to the value 
corresponding to a high-fitness model of the inter- 
action. 
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Figure [7] reveals a qualitative difference in the 
way the specified parameters affect the motion of 
the Magellanic Clouds. The influence of the LSR 
circular velocity is significantly stronger than the 
role of the flattening of the Galactic potential. 
This is not surprising. While the flattening pa- 
rameter only has an impact on the distribution of 
mass of the Galaxy, variations in the LSR circu- 
lar velocity also lead to changes of the LMC/SMC 
galactocentric velocities. Moreover, these effects 
amplify each other regarding their impact on the 
orbits of the Clouds (see Section [5|). If 60 is in- 
creased, the total mass of the Milky Way increases 
implying a higher gravitational potential energy of 
the Clouds with respect to the Galaxy. Regarding 
Equation (fl7|) and the related discussion, the in- 
creasing LSR circular velocity reduces the galac- 
tocentric velocities of the Clouds, i.e. their ki- 
netic energy. Thus, the total energy of the LMC- 
Galaxy and the SMC-Galaxy pairs changes with 
0o rapidly. 

The lower row of Figure [7] shows that the or- 
bits of the Magellanic Clouds projected to the 
plane of sky are always offset from the position 
of the Magellanic Stream, regardless of the choice 
for the [q, G>o] pair. Closer look at Figure [15] re- 
veals that such a deviation is controlled by the 
northern component fj}™ of the LMC proper mo- 
tio n and cannot be remove d unles s the results 
by Kallivavalil et al.l (|2006bl laT) and Piatek et al 



(2008) are revisited. The particles forming the 
leading and trailing streams in our restricted N- 
body model (and generally in every tidal model) 
essentially follow the orbits of their progenitors. 
This purely tidal approach is thus limited. As long 
as the LMC/SMC projected orbits are offset from 
the observed Magellanic Stream, it is difficult to 
model the full extent of the large-scale structure 
of the Stream in the position-position-LSR radial 
velocity space (see Section 16. 8j) . 

IChoi et al] ( 2007 ) have demonstrated that tidal 
tails originating in massive satellite galaxies be- 
come offset from their progenitors' orbits to some 
level. However, unless additional physical pro- 
cesses introduce the dissipation of energy and an- 
gular momentum of the particles in the tidal tails, 
such a phase-space deviation between the stream 
and the orbit of its progenitor is limited. Figure [8] 
shows how this is reflected by the kinematics of 
the Magellanic Clouds and the Stream. The high- 



fitness models are split into the groups of / > 0.60 
and / < 0.50, respectively. It is a common feature 
of the best models (/ > 0.60) that the LSR ra- 
dial velocity along the LMC orbital track agrees 
well with the famo us linear velocity p rofile of the 
Magellanic Stream (|Bruns et alj|2005| ). 

Despite the previously discussed behavior of 
tidal models, the apparent offset of the recent 
(-300 Myr < t < 0) LMC/SMC orbits from the 
Magellanic Stream was reproduced partially by 
our simulations. It did not occur due to a large 
deviation between the tidal particle tail and the 
orbits of the Clouds. The modeled projection of 
the Magellanic Stream to the plane of sky is still 
well aligned with the previous revolutions of the 
LMC/SMC orbits (see Figure |TS|) . but the preced- 
ing orbital cycles are projected to different posi- 
tions in the plane of sky reflecting the geometry 
and physical properties of the LMC-SMC-Milky 
Way system. 

The subsequent revolutions of the orbits of the 
Magellanic Clouds are shifted with respect to each 
other in our good models of the interaction with 
the Milky Way. This is to be attributed to the 
projection effect of the configuration of the system 
involving the observer (Sun), the Galactic Cen- 
ter and the Magellanic Clouds (Figure [2} . The 
flattening of the Galactic potential must be taken 
into account, as it prevents the LMC from hav- 
ing an orbit confined within a 2D plane. Figure [5] 
shows the intensity of the radial and the axial com- 
ponents of the gravitational field of the logarith- 
mic halo. The ratio of the components depends 
strongly on the position [R, z] with respect to the 
Galactic Center. That leads to the formation of 
a 3D orbit of the LMC. We did not mention the 
SMC, because its low mass makes it dependent not 
only on the gravitational field of the Galaxy, but 
also on the attraction by the LMC. 

The previous studies of the dynamical evolu- 
tion of the Magellanic Clouds assumed a good 
alignment of the LCM/SMC orbits with the Mag- 
ellanic Stream in the p osition-position-LSR radial 
velocity space (see, e.g. Gardiner fc Noguchilll996t 



Mastropietro et al.l 120051 IConnors et al.l 12006) 



However, such an assumption is at odds with 
the HST proper motion measurements for the 



Clouds (Besl a et al.l 12007). We have seen that 



nothing chang es even for the revi s ited v alues of 
6 and fl Q bv iReid fc Brunthalerl (|2004h . What 
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it represents to modeling the large-scale features 
of the Magellanic Clouds was briefly discussed in 
this section and it will be further addressed later 
in Sections 16.81 and 16.91 

6.5. LMC-SMC encounters 

The Magellanic Clouds sha re a common low 



density gaseous envelope (see iBriins et al.l |2005) 



which is usually considered one of the signs in- 
dicating that the Clouds have formed a grav- 
itationally bound couple even for several Gyr. 
Gravitational binding of the Clouds lasting sev- 
eral Gyr turned out to be extremely rare in our 
si mulations. However, the proper motion data 
by Kallivavalil et al. ( 2006blla ) and bv lPiatek et al. 
(|2008l ) imply that the distance between the LMC 
and the SMC is very likely to have been shorter 
than today for several hundred Myr. This may 
be sufficient to explain the presence of the HI 
envelope, as our simulations have shown that an 
envelope surrounding the Magellanic Clouds may 
be formed rather quickly on the specified time- 
scale. 

The information about the LMC-SMC distance 
in the recent past comes directly from the combi- 
nation of Figures [10] and [TTJ We have calculated 
the scalar product of the present galactocentric ve- 
locity vectors of the Magellanic Clouds for all the 
high-fitness models. Figure [TU] shows the fitness 
of these models as a function of the actual value 
of the scalar product. The current LMC and SMC 
galactocentric velocity vectors are nearly parallel 
to each other in all cases and their angular devi- 
ation does not exceed 20° with the mean value of 
only sa 10°. It is natural to ask whether the spa- 
tial separation of the Clouds is increasing or de- 
creasing at present. Figure [TT] offers the answer. 
If the time dependence of the distance between 
the Clouds is plotted over the last 700 Myr for the 
models of / > 0.55, one may clearly see that the 
LMC and the SMC have been receding from one 
another for rj 100 Myr. 

To quantify the rate at which the LMC-SMC 
separation Z?i_ s changes we have calculated the 
function dD\^ a /dt at the present time t = t$ = 
Gyr. If the function dD\_ s /dt is evaluated for 
the models plotted in Figure [TTJ it returns the 
values between 50kpcGyr _1 and 120kpcGyr _1 . 

The preceding paragraphs lead to the conclu- 



sion that the orbital history of the Magellanic 
Clouds in the high-fitness models always involves 
a close encounter between the LMC and the SMC 
at the time of -250 Myr < t < -80 Myr. It 
is quite interesting to see whether such behavior 
is outstanding or typical within the LMC/SMC 
proper motion ranges established by the recent 
HST measurements. 

Figure [T^] illustrates how the rate at which 
the current spatial separation of the Magellanic 
Clouds changes depends on the proper motion of 
the SMC. We have examined two cases. In one 
case, the proper motion components of the LMC 
were fixed to the values of a selected very good 
model (/ > 0.60). The model was chosen in such 
a way that its LMC proper motion components are 
close to the midpoint values of the intervals estab- 
lished by Equations ([T]). The SMC proper motion 
was varied within the entire ranges for /i™ c and 
^J^ c involved in our study. In all cases, the Mag- 
ellanic Clouds were found to be receding from one 
another (Figure [T2l left plot), i.e. 



df 



> ( smc smc 



)>o. 



(20) 



In the other case, we have kept the values of all the 
parameters equal to those ones of the mentioned 
high-fitness model, but the LMC proper motion 
was modified to reach the minimum galactocen- 
tric velocity. Calculating the rate dD\_ s /dt then 
yielded the result depicted in the right hand plot 
of Figure 1121 Obviously, no change occurred and 
Equation (gDJ) still holds. 

Thus, we have seen that strong support exists 
for a remarkable feature of the kinematics of the 
Magellanic Clouds. The recent estimates of the 



LMC/S MC proper motions by jKallivavalil et al. 
(|2006bl fah and by IPiatek et al.l (|2008l ) are very 



likely to introduce a close (A-s ~ lOkpc) past 
encounter between the Magellanic Clouds at the 
time of ps —100 Myr for an arbitrary choice for 
the LMC/SMC proper motions. 

It is reasonable to expect a long-term gravi- 
tational binding between the Clouds of the order 
of 10 9 yr to be a natural condition for repeated 
encounters in the LMC-SMC system. The corre- 
sponding alignment of the past orbital tracks of 
the Clouds in the phase space requires a very spe- 
cial setup of the parameters for the model of the 
interacting system, making such a case very un- 
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likely. However, Figure [13] reveals a surprising pic- 
ture of the LMC/SMC orbital motion in the deep 
past. If the separation between the Clouds for the 
best models (/ > 0.60) is plotted as a function of 
time, similar behavior is found in all cases. After 
the previously discussed close approach a couple 
of Myr ago, the mutual distance of the Magel- 
lanic Clouds reaches 150 kpc, exceeding the value 
of 200 kpc for some configurations. The Clouds 
become unbound shortly after their most recent 
encounter. Nevertheless, every model contains a 
second close approach at the time t < —2.5 Gyr 
with the separation of the Magellanic Clouds drop- 
ping below « 15 kpc. As we will show later, this 
encounter is the event triggering the formation 
of an extended trailing tidal tail (the Magellanic 
Stream) and of its leading counterpart (the Lead- 
ing Arm). 

6.6. Redistribution of mass in the Magel- 
lanic System 



We have stated in Section I6TT1 that only a sub- 
set of the models identified by the genetic al- 
gorithm in the parameter space may be consid- 
ered satisfac tory reproduction s of the observa- 



tional data ( Briins et al. 2005). In general, an 



acceptable model of the interaction between the 
Magellanic Clouds and the Milky Way is supposed 
to reproduce the extent of the observed Magellanic 
Stream in the 3 D space composed of the LSR ra- 
dial velocity and of the position in the plane of 
sky. The second large-scale HI structure associ- 
ated with the Clouds - the Leading Arm - is not 
used to separate good models from the unsuccess- 
ful ones. 

The Leading Arm appears in tidal models of 
the LMC-SMC-Galaxy interaction as a natu- 
ral counterpart of t he Magellanic Stream (se e 
Gardiner fc Noguchil Il996t IConnors et ahl l2006h . 
but its acceptable reproduction by the means of 
numerical simulations remains a challenge. Our 
parameter study has not changed that picture - 
the modeled distribution and kinematics of HI 
in the Leading Arm region and around the main 
LMC and SMC bodies remains similarly unsatis- 
factory over the parameter space, especially fail- 
ing to reproduce the observed morphology of the 
Leading Arm. In terms of the fitness function, the 
value of / never exceeds ~ 0.35 if calculated only 
for the Leading Arm. 



Reliable modeling of the evolution of the cen- 
tral regions of the Magellanic Clouds currently 
eludes the level of sophistication of our numerical 
code. The inner parts of the Clouds contain high- 
density baryonic mass and in order to study their 
evolution a detailed model involving self-gravity 
and hydrodynamical processes which account for 
the mass exchange cycle between stars and gas, 
is required. Our restricted N-body code overesti- 
mates the amount of gas tidally stripped from the 
LMC/SMC central regions because the tides are 
not balanced by self-gravity and the dissipative 
behavior of gas. As a consequence, the modeled 
column density is lower than observed values in the 
inner 5 kpc of the Clouds. This effect is stronger in 
the case of the SMC as it is the dominant source of 
matter for the Magellanic Stream and the Leading 
Arm (see Figures fT4l or [T6|) . The particle distribu- 
tion in the LMC is affected by tidal heating by the 
Milky Way and the original particle disk is turned 
into a 3 D spherical structure. 

6.7. Definition of a good model 

The Magellanic Stream turned out to be the 
appropriate indicator of the quality of our mod- 
els. Its formation is very sensitive to the choice of 
the model parameters and critically influences the 
resulting fitness. Figure [T4l illustrates these facts. 
While the model of / = 0.61 is able to fit the 
basic features of the Magellanic Stream both in 
the projected H I distribution and the LSR radial 
velocity profile, a typical simulation representing 
the model group of / < 0.50 places the Magel- 
lanic Stream to the position-position-LSR radial 
velocity space incorrectly. The Stream is extended 
insufficiently both in position and the LSR radial 
velocity (upper row of Figure [14]), and slope of the 
simulated profile of the LSR radial velocity along 
the Stream exceeds the observed one already at 
the Magellanic longitude of ss —30° (lower right 
plot of Figure IT4]) . Generally speaking, the de- 
scribed behavior of the modeled trailing stream 
is responsible for the resulting fitness of a given 
model, and was used to define the desired thresh- 
old level of the fitness value. 

The obvious lack of models within the fitness 
range of (0.50,0.55), clearly visible in Figures [4] 
and [TO] was combined with the specified features 
of the high-fitness models. This yielded the de- 
sired threshold level of the fitness function quite 
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naturally. Figure [4] shows that the fitness of the 
model establishing the upper limit to the men- 
tioned gap in the distribution of the fitness val- 
ues lies slightly below 0.55 (it is 0.546 actually). 
However, the transition between the good and un- 
satisfactory models is always gradual, as we have 
seen, and thus the value / = 0.55 can be taken for 
the threshold value of the fitness function with no 
loss of generality. 

6.8. Morphology and kinematics of the 
Magellanic Stream 

In order to discuss our results regarding the 
reproduction of the large-scale distribution of 
HI associated with the Magellanic Clouds, we 
compared a selected high-fidelity model with the 
lo w-resolution d ata-c ube of the HI observations 
bv I Brims et al. ( 2005 ). The simulated 3D data- 
cube depicted in Figure [TH] (already shown in Fig- 
ure IT4"]) was produced by a model with the fol- 
lowing values of the LSR circular velocity and the 
flattening of the Galactic halo: = 232 km s -1 , 
q = 0.81. These values yield the total enclosed 
mass of the Milky Way within the radius of 250 kpc 
of 2.20- 10 12 M©. 

Assuming the recent measurements of the LMC 
and the SMC proper motions, our tidal model was 
able to produce a trailing tail similar to the ob- 
served Magellanic Stream regarding both the LSR 
radial velocity profile and the extent in the pro- 
jected position. While the LSR radial velocities 
measured along the Magellanic Stream are repro- 
duced very well, the far tip of the simulated stream 
at the position of b > —60° and I ps 70° is dis- 
placed compared to the observations (Figure [T5J . 
Such a displacement appears in every high-fitness 
model identified by the genetic algorithm in the 
studied parameter space. 

The column density of HI in the Magellanic 
Stream decreases rather smoothly towards its far 
tip. This feature exists also in our high-fidelity 
simulations of the tidally induced formation of 
the large-scale Magellanic structures. The mid- 
dle row of Figure [15] shows how the mean HI 
column density on the line perpendicular to the 
gr eat circle of the Magella nic longitude (defined 
in lWannier fc Wrixon 1972 ) changes with the lon- 
gitude. The simulated profile of the column den- 
sity (right hand plot) has a steeper slope than the 
observed Stream and its gradual decrease ends at 



the Magellanic longitude of —70° compared to 
~ —90° for the observations. 

The plots in the lower row of Figure [15] also of- 
fer the view of the area in the plane of sky where 
the H I structure ca l led th e Leading Arm is lo- 
cated. iBriins et al.l ( 2005 ) introduced the divi- 



sion of the structure into three systems of gaseous 
clumps named LAI, LA II and LA III. Their ap- 
proximate positions in the Galactic coordinates 
[b,l] are [-15°, 300°] (LAI), [15°, 290°] (LA II), 
and [15°, 270°] (LA III). The simulation shown in 
Figure [15] was able to reproduce the part LA III 
and a fraction of the clump LA II. Some of our 
high-fitness models offer quite a good match of 
the overall kinematics and spatial extent for both 
observed structures LA II and LA III. Neverthe- 
less, no model reproducing both the system LAI 
and the Magellanic Stream at the same time was 
found. 

6.9. Formation of the Magellanic Stream 

It is interesting to take a look at the compo- 
sition of the simulated Magellanic Stream, and at 
the process of its build-up over the course of time. 
The lower left plot of Figure Q3] reveals that most 
of the particles forming the Magellanic Stream in 
the selected high-fitness model are of the SMC 
origin. A number of the LMC particles are spread 
over the full extent of the Stream as well. This is 
a typical result and it is based on the fact that the 
LMC particles are more tightly bound due to the 
significantly higher mass of the LMC compared to 
the SMC - the LMC:SMC mass ratio is 10:1 in 
our simulations. 



It h as already been discussed by IConnors et al 



(2006) ;hat the Magellanic Stream might consist 
of several filaments projected to the same region 
in the plane of sky. Such filaments would be the 
tidal debris spread along the past orbits of the 
Clouds. Our study has confirmed such a scenario. 
Figure [15] shows the projected orbits of the Mag- 
ellanic Clouds from t = — 3Gyr to the present, 
plotted over the contours of the integrated column 
density of HI (lower right plot). The SMC has 
crossed the area currently occupied by the Magel- 
lanic Stream three times, while two passages have 
occurred in case of the LMC. If we disentangle the 
projected orbital paths, we will find out that the 
previous orbital revolutions appear more diagonal 
in our plot. They are placed to lower Galactic lon- 
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gitudes. The resulting trailing tail is thus a mix- 
ture of particles released from the LMC/SMC due 
to the tidal stripping during several epochs cov- 
ering the time range from present to w — 3.3Gyr 
when the redistribution matter was triggered by a 
close LMC-SMC encounter. 

The complex structure of the modeled Stream 
can be seen in the lower plot of Figure Ql)] The 
LMC/SMC particles are color coded according 
to the epoch when they were stripped from the 
Clouds. The Stream and the Magellanic Bridge 
are dominated by the SMC particles, while most 
of the LMC particles remain in the vicinity of the 
main LMC body and have been gradually tidally 
heated due to the gravitational field of the Milky 
Way. 

The past orbital history of the Magellanic 
Clouds for the discussed high-fitness model can 
be seen in the lower plot of Figure [16] For the 
last 4 Gyr the Clouds have moved within the 
Galactocentric radius of 200 kpc, undergoing two 
close (D\_ s s» 10 kpc) encounters at the times of 
—3.3 Gyr and —0.15 Gyr, respectively. The im- 
pact of the encounters on the distribution of the 
LMC/SMC particles was different. 

To quantify the influence of the encounters in 
the LMC-SMC-Galaxy system, we have counted 
the number of particles shifted from the original 
orbits of the radius Ti to the distance of 2r^ from 
the center of their progenitor. If this function 
is evaluated for the time bins of 500 Myr from 
t = —4 Gyr, we obtain the upper plot of Fig- 
ure HU Apparently, the number of the disturbed 
particles was increased already in the time bin 
of (—3.5,-3.0) Gyr as a consequence of the first 
LMC-SMC close approach. This level was later 
only mildly decreasing, as more of the disturbed 
particles were stretched in the tidal field of the 
Milky Way, satisfying the above stated definition 
of strongly affected particles. 

The second close LMC-SMC encounter was sig- 
nificantly more dramatic, disturbing more than 
10% of all particles. This event is the origin of the 
common envelope surrounding both the Clouds, 
and of the turbulent gaseous filament connecting 
the main bo dies of the Clouds, named the Magel- 
lanic Bridge (| Brims et al. 1 12005ft . Substantial frac- 
tion of the particles contributed to the Magellanic 
Stream as well. 



The results summarized in the preceding para- 
graphs and depicted in Figure [16] have a remark- 
able consequence regarding the formation history 
of the Magellanic Stream. Over the last ~ 3 Gyr 
it was supported at an almost constant rate by 
the gas tidally stripped from the Clouds. How- 
ever, this rate was increased by a factor ps 20 
due to a recent encounter between the Clouds at 
t = —0.15 Gyr. We would like to mention that on 
the qualitative level these results are an appropri- 
ate description of all the high-fitness models found 
in the parameter space for the LMC-SMC-Milky 
Way interaction. 

We have found that the high-fitness simulations 
prefer the oblate (q < 1) shape of the logarithmic 
halo of the Galaxy over the prolate configuration 
(see Figured]), but several reasonable models of q 
as high as 1.10 exist. As long as good models of 
the interacting system exist for the prolate halo, 
this shape of the Galactic gravitational potential 
cannot be ruled out. However, we have shown that 
such a possibility is substantially less likely than 
the case of q < 1. Moreover, only two acceptable 
models with significantly prolate (q > 1.10) DM 
haloes of the Milky Way have been identified for 
a reasonable total mass of the halo, i.e. below the 
limit of « 3 • 12 M Q (see Section l3~Tj) . 

The impact of the flattening parameter on the 
phase-space positioning of the LMC/SMC orbits 
has been addressed in Section 16.41 of this paper. 
The shape of the Galactic halo also influences the 
redistribution of the matter associated with the 
Magellanic Clouds due to the resulting tidal field 
of the Milky Way. This effect is the key to the pref- 
erence for the oblate configuration of the DM halo 
of our Galaxy. We have calculated the tidal force 
exerted by the Milky Way on two points separated 
by Ar = 10 kpc and moving around the Galactic 
Center on circular polar orbits of the radii 50 kpc, 
100 kpc and 150 kpc. Figure [IT] shows the result. 
In general, the oblate configuration offers the most 
efficient tidal stripping, improving the probability 
that a sufficient redistribution of the LMC/SMC 
particles occurs within the available time-scale of 
several Gyr. 

7. Summary and Conclusions 

We have performed a detailed exploration of 
the parameter space for the interaction involv- 



15 



ing the Magellanic Clouds and the Milky Way. 
T he method app l ied w as similar to the approach 
bv lRuzicka et al.1 (|2009f ). but the results were quite 
different as a new free parameter was involved. 

The recent work by IShattow fe Loebl <|2009h 
lead to the conclusion that the modified values 
of the LSR circular velocity and of the LSR an- 
gular rotation rate (jReid fc Brunthalerl 120041 ) in- 
troduce a significant change to the orbital his- 
tory of the Magellanic Clouds. We have followed 
their study in order to show whether the large- 
scale distribution of HI in the Magellanic sys- 
tem can be simulated successfully even for the in- 
creased values of the LMC/SMC prope r motio ns 
measur ed recently b y Kallivavalil et al. ( 2006bl la) 
and bv lPiatek et al.l ( 20081) if the revisited view of 
the galactic rotation (jReid fc Brunthalerl 120041 ) is 
adopted. 

The study by iRuzicka et all (|2009h ruled out 
the tidal models of the LMC-SMC-Galaxy inter- 
action because they were unable to reproduce the 
basic kinematic and morphological features of the 
most prominent structures originating in the Mag- 
ellanic Clouds, i.e. of the Magellani c Stream and 
the Leading Arm ( Brims et al. 20051 ) . Their auto- 
mated analysis of the parameter space for the in- 
teraction has shown clearly that the recent proper 
motion measurements imply the past LMC/SMC 
orbits that allow for efficient tidal stripping of gas 
from the Clouds. 

Unlike Ruzicka et al. (2009). we have taken into 
accou nt the observational data bv lReid fc Brunthaler 
(|2004l ) and treated the LSR circular velocity as 
a free parameter varied withi n the range 0n = 
(210, 2 60) kins' 1 . Following iReid fc Brunthaler 
(2004), a modified value of the LSR angular ro- 
tation rate fio of 29.45 kms -1 kpc -1 was adopted 
as well, yielding the studied range of the Solar 
galactocentric distance Rq of (7.13, 8.83) kpc. As 
shown in Section [5] of this paper, the new values 
of 9o and f^o affect the galactocentric positions 
and velocities of the Magellanic Clouds, as well as 
the mass of the Galaxy. 

The exploration of the modified parameter 
space for the interaction has revealed a remark- 
able qualitative change regarding the features of 
the resulting candidates f or acceptable tida l mod- 
els. The expectations of Shattow fc Loebl (2009) 
have been confirmed because satisfactory repro- 
duction of the large-scale Magellanic structures 



was possible for the LMC/SMC proper motion 
data by the HST. Moreover, such quality mod- 
els have been localized for a large number of the 
proper motion combinations (Figure 2|) . This is 
due to the fact that the proper motions no longer 
play the exclusive role in establishing the actual 
3 D galactocentric motion of the Clouds. 

The LMC/SMC velocity (and position) vectors 
are determined not only by their proper motions 
but also by the LSR circular velocity and the LSR 
angular rotation rate. The current phase space co- 
ordinates of the Clouds are linked with the mass 
distribution of the Galaxy due to 0o and Qq. As a 
consequence, a great variability and freedom have 
been introduced in the parameter space concern- 
ing the options to choose the orbits of the Magel- 
lanic Clouds allowing for the efficient tidal redis- 
tribution of mass resulting in the formation of the 
Magellanic Stream. 

All the high-fitness models of the interac- 
tion localized in the studied parameter space in- 
volved two close LMC-SMC encounters within 
the last 4 Gyr. The first one occurred at the time 
t < —2.5 Gyr and triggered the evolution of the 
Magellanic Stream. This encounter caused the 
tidal heating of the outer regions of the origi- 
nal LMC/SMC particle discs. Subsequently, the 
disturbed particles were spread along the orbital 
paths of the Clouds due to the tidal stripping by 
the gravitational field of the Galaxy. 

The latter encounter was placed only as re- 
cently as -150Myr (-250 Myr < t < -80Myr), 
but its impact on the LMC/SMC particles was 
at least by a factor of 10 stronger compared to 
the first encounter (for more rigorous discussion 
see Section This event was also the begin- 

ning of the formation of the filament connecting 
the Clouds, composed of turbulent g as and young 



stars , called the Magellanic Bridge (| Brims et al 



2005). It is a notable fact that the recent en- 



counter between the Clouds is a rather general 
feature of the interaction. Figure [12] demon- 
strates that it is likely an intrinsic property of 
th e recent proper motion data. The studies 
by_ Kallivavalil et al.1 ( 2006bl fah and bv lPiatek et al 



(|2008l ) seem to introduce a LMC-SMC approach 
to the distance of w 10 kpc at the time of 
« -150 Myr. 
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8. Discussion 

Although it is usually assumed that a long- 
term gravitational binding between the Clouds has 
exi sted, we have confirmed th e previous findings 
by iRuzicka et all (|2007l l2009h that such a condi- 
tion is not necessary regarding the tidally induced 
redistribution of mass associated with the Clouds 
leading to the formation the Magellanic Stream. 
Figure Q2] shows clearly that after the first en- 
counter, the spatial separation of the Clouds usu- 
ally exceeds 200 kpc. 

The assumption of the Clouds being gravita- 
tionally bound is supported by several more obser- 
vational indications. First of all, the Magellanic 
Clouds are surrounded by a low-density gaseous 
envelope which is a natural feature associated with 
two interacting bodies sharing a common history. 
However, we have found that the second of the 
mentioned LMC-SMC encounters is able to pro- 
duce such a diffuse gaseous structure. 

The gravitational binding of the Clouds is also 
often substantiated by their unique composition 
which is quite outstanding within the neighbor- 
hood of the Milky Way. Most of the satellites of 
the Galaxy are dwarf spheroidal while the LMC 
and the SMC are highly unevolved gas-rich irregu- 
lar galaxies. The likelihood that such a couple was 
formed by chance in the Local Group of Galaxies 
is very small, obviously. These issues are serious, 
but they cannot be addressed within the frame- 
work of our parametric study. 

It is one of the most surprising results of the ge- 
netic algorithm-based exploration of the parame- 
ter space that good models have been localized 
independently of the total mass of the Milky Way 
(see Figure [5]). This behavior of the interacting 
system occurred due to the LSR circular velocity 
that links the present phase-space positions of the 
Clouds with the actual enclosed mass of the galac- 
tic DM halo. The consequence of such a complex 
nature of the system is nicely manifested by Fig- 
ure [6] 

However, even the remarkable reduction of 
the uncertainties of the LMC /SMC prop er mo- 
tio ns achieved b y iKallivavalil et al. | (12006 b, and 
by IPiatek et al.l ( 2008 ) has left error boxes large 
enough to accommodate virtually every request 
for the proper motions to compensate the actual 
choice for the LMC circular velocity. Hence, the 



current phase-space positions of the Clouds can 
be selected in order to make the tidal model work 
over a wide range of the enclosed Galactic mass. 
Reasonable simulations of the Magellanic Stream 
have been made over the range of 0.6 • 10 12 M Q to 
3.0 • 10 12 M for the total mass of the Milky Way 
enclosed within the radius of 250 kpc. 

It cannot be omitted that as the Galactic mass 
decreases due to the decrease of the LSR circu- 
lar velocity, the density profile of the DM halo 
steepens in the inner region of the Galaxy. This is 
caused by the rescaling of the rotation curve of the 
Milky Way yielded by the relation J7o = Qo/Ro = 
const. Such behavior is helpful in respect to our 
goals, but it must be treated with care. 



Besla et al.l (|2007l) have pointed out that vari- 



ous observations have put constraints on the upper 
limit for the total mass enclosed within w 50 kpc 
from the Galactic Center. However, taking such 
findings into account does not disprove our results. 
The limits on the total mass in the inner part of 
the Milky Way can be considered simply as addi- 
tional constraints reducing the number of accept- 
able parameter combinations. 

The redistribution of the matter associated 
with the Magellanic Clouds followed a sc enario 
similar to that one by Gardiner et al.l ( 1994 ) . The 
formation of the Magellanic Stream (and of its nat- 
ural counterpart - the Leading Arm) was triggered 
by a LMC-SMC encounter 2.5 Gyr ago. However, 
the Stream as old as sa 2 Gyr is at odds with the 
rec ent conclusions by 
bv iBland-Hawthorn ( 



Stanimirovic et al. I (120081) or 
2009). They argue that the 



interaction of the relatively cold H I of the Stream 
with the hot ambient halo of th e Milky Way results 
in th e thermal fragmenta tion (jStanimirovic et all 
20081 ) and the ablation (jBland-Hawthornl |2009h 



of the HI clouds. Such a decay yields the maxi- 
mum survival time of the H I clouds not exceeding 
~ 1 Gyr. However, our model introduces a con- 
tinuous replenishment of the Stream gas, and a 
strong boost to this process due to the second 
LMC-SMC encounter at « -0.15 Gyr. 

We also have to point out the fundamental 
difference regarding the phase-space structure of 
the Magellanic Stream if the tidal and ram pres- 



sure m odels are compared. iMastropietro et al 
([20051 ) have shown that the ram pressure-induced 
Stream possesses a rather compact phase-space 
structure. On the other hand, we have seen 
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that the tidal mechanism leads to the Magellanic 
Stream composed of several filaments of different 
ages. These filaments are the remnants of the 
past LMC/SMC orbital revolutions and might lie 
in significantly different distances even as large as 
15 kpc (see Figures [TBI and IT6| . Since the study 
bv IStanimirovic et al.l (|2008h assumed the Stream 
to lie not further than 50 kpc from the Galactic 
Center, their results cannot fully be applied to 
our models. The decay of the HI clouds due to 
the interaction with the hot gaseous halo of the 
Milky Way might be slower at larger radii as the 
density of the ambient Galactic medium falls with 
the galactocentric distance. 

It is a challenging problem for the tidal models 
of the Magellanic System to explain the absence 
of stars in the Stream. Whilst evolutionary sce- 
narios taking ram pressure stripping as the dom- 
inant process in the formation of the Magellanic 
Stream explain the problem naturally, tides affect 
both stars and gaseous clumps. Hence, the miss- 
ing stars pose a serious issue regarding the tidal 
origin of the Stream. 

The study bv lBekki fc Chibal (|2009l ) has shown 
by the means of self-consistent numerical simu- 
lations that the key parameter determining the 
stellar content in a tidal model of the Magel- 
lanic Stream is the relative extent of gas and 
stars in the SMC. Assuming that the radii of the 
HI and stellar distr ibutions in the SMC are of 
2 < r H i/r sta rs < 4, iBekki fc Chibal (|2009t ) have 
successfully reproduced the observed composition 
of the Magellanic Stream. The corresponding 
stellar distribution was disturbed by the recent 
(f « -200 Myr) LMC-SMC encounter and con- 
tributed to the extended SMC/LMC halos, but 
no stars appeared in the Magellanic Stream. The 
assumption on rHi/r s t a rs has a solid observational 
support, as the HI diameters of gas-rich galaxies 
are observed to b e significantly larger than the ir 
optical disks (e.g. Broeils fc van Woerden 1994) . 

Stars and gaseous clumps are treated identically 
as test particles in our model and they satisfy the 
same equations of motion. Whether a given parti- 
cle represents a star or an H I cloud depends on its 
initial distance from the SMC (LMC) center and 
on the assumed radius of the stellar body. This 
quantity was estimated through the use of obser- 
vational and theoretical constraints summarized in 
the previous paragraph and was set to 0.5rdi s k (see 



Table [J). 

The recent LMC-SMC encounter strongly af- 
fected the test particles w ithin the radius of the 
stellar body, but similar to Bekki fc Chib a (2009) 
most of the particles contributed to the common 
envelope of the Clouds. However, < 5 % of such 
particles were moved to the Stream and may thus 
be considered its stellar contamination. Unfortu- 
nately, this is inevitable in a restricted N-body 
simulation. Our model overestimates the amount 
of matter stripped from the SMC center due to 
its over-simplified description, where self-gravity 
and dissipative gas dynamics will act against the 
tidal stripping as significant restoring forces (see 
Section 16.61) . 

9. Open questions 



Th e values of the LSR circular velocity bv lReid fc Brunthaler 
certainly represent a progress towards the 
resolution of the difficulties arising from the recent 
HST proper motions of the Magellanic Clouds. 
However, we have not revealed a parameter set 
allowing for the correct modeling of the far tip of 
the Magellanic Stream. An offset of a model from 
the observations was present (Figure 1X5"]) . This 
occurred because the LMC/SMC orbits have not 
crossed the corresponding region of the position- 
position-LSR radial velocity space. Nevertheless, 
we should keep in mind the simplicity of our 
model. Neither the dissipative hydrodynamical 
processes were included, nor did we account for 
th e alternative s c enario s such as the recent idea 
bv lNidever etaD (|2008h . 

They have analyzed the kinematics of the Mag- 
ellanic Stream and the Leading Arm, with the fo- 
cus on the transition regions between these struc- 
tures and the main bodies of the Clouds. Their re- 
sults suggest that bo th the Magellanic S tream and 
the structure LAI ( Brims et al. 20051 ) might be 
evolutionarily related to the region in the south- 
east of the LMC where a massive star formation 
takes place. It is possible that such a star forming 
activity has caused a strong blowout of gas which 
was later redistributed by the tidal or ram pres- 
sure stripping. 

The scenario proposed bv lNidever et al. (2008) 
might be the complementary process to the tidal 
stripping concerning the origin of the Leading 
Arm. While the tidal model succeeded reproduc- 
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ing the clumps LA II, LA III, the gas blow-out 
would transport gas to the position of the complex 
LAI. In general, such a process is able to eject gas 
with an arbitrary direction of its momentum with 
respect to the LMC motion. It might allow for fill- 
ing the regions of the phase-space unreachable by 
the tidal/ram pressure models with gaseous mat- 
ter. 
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Fitness function 



We have further improved the fitness function devise d by Ruzicka et al.1 (2007 ) in order to handle the 
weaknesses of their scheme revealed later in the work bvlRuzicka et all (12009^ The first component of our 
fitness function is still based on the approach by iTheis &: Kohlei (|200lh who proposed a generally applicable 
technique of comparing the relative intensities of the corresponding pixels in the modeled and observed 
data-cubes. 

Both modeled and observed HI column density values are scaled relative to their maxima to introduce 
dimensionless quantities. Then, we get 



f 



N„ ■ N x ■ N, 



N v N y N x 

EEE 

y <=i j=i k=i 



1 



obs 

ijk 



a 



mod 
ijk 



(Al) 



where cr° b fe s , cr^ d are normalized column densities measured at the position [j, k] of the i-th velocity channel 
of the observed and modeled data, respectively. N v is the number of separate LSR radial velocity channels 
in our data. (N x ■ N y ) is the total number of positions on the plane of sky for which observed and modeled 
H I column density values are available. 

The second component of the fitness function performs a search for structures in the position-position- 
velocity space of the data-cube. The modeled distribution of HI is compared to its observed counterpart 
regardless of the exact levels of the H I radio emis sion. It combines the enhancement of structures in the data 
by their Fourier filtering (see lRuzicka et alj|2007|) with the subsequent check for empty/non-empty pixels in 
both data-cubes. The corresponding component of the fitness function is defined as follows: 
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where pix° b £ £ {0, 1} and pix™£ d € {0, 1} indicate whether there is matter detected at the position [i,j, k] 
of the 3 D data-cube of the observed and modeled Magellanic System, respectively. 



Ruzicka et al.l ([20071 ) showed that the search for structures significantly improves the performance of the 



genetic algorithm if the structures of interest occupy only a small fraction of the system's entire data-cube 
(< 10% in the case of the Magellanic Stream and the Leading Arm). 

The fitness function /2 is a measure for the agreement of the shape in the data. No attention is paid to the 
actual H I column density values. However, the efficiency of such a search for structures becomes lower with 
the increasing resolution of the compared data-cubes, as the relative number of non-empty pixels usually 
decreases in such a case, and the function fa returns very small values. To resolve the described difficulty, 
we have proposed a new component of the fitness function which is based on Equation (|A2[) . In fact, we 
have introduced another simplification level for the view of the properties of the studied system. Similarly 
to Equation (|A2[) the number the overlapping non-empty pixels in the observed and modeled data-cubes is 
counted. This procedure is performed repeatedly for a series of N data-cubes of an increasing resolution. It 
is significantly easier to achieve a good match in case of low-resolution data, but the sensitivity to structures 
is rather poor, and only global positioning of the system as a whole is evaluated. High-resolution data 
provide more detailed information on the system, but the corresponding complexity prevents the comparison 
from driving the genetic algorithm efficiently. The newly designed component of the fitness function may be 
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described as 



JV 

E 

m=0 



JV™ N™ 

EEI 

i=l j=l k=l 



E 2 "'EE Ep^f-pi^ 



A? W? 1 IV" 



(A3) 



E2-EEE i'i<r 



t=l J=l fe=l 

where iV™, -/V™, and AT™ are the dimensions of the m-th data-cube and 



N™ < AC" +1 A A/" 1 < N" l+1 A AC™ < N. 



m+l 



Obviously, we have introduced a weighting factor of 2 m in Equation (|A3I) . It reflects the above discussion 
and strongly emphasizes the models that are able to satisfy the high-resolution observational data. 

Ruzicka et al.l (|2007l ) also recommended and successfully applied a system-specific comparison. In the 



case of the Magellanic Clouds, the very typical linear radial velocity profile of the Stream including its high 
negative velocity tip was considered important. The slope of the LSR radial velocity function is a very 
specific feature, strongly dependent on the features of the orbital m otion of the Cl o uds. We have slightly 
modified the original definition for the LSR radial velocity check by iRuzicka et al. (2007), and the fourth 
fitness function component is defined as 



h = 



1 



E 



pix? bs (n m i n ) — pix™° (u m i n ) 



(A4) 



where pix° bs (u 



and pix™ od (u m i n ) are the pixels with the minima of the observed LSR radial velocity 
profile of the Magellanic Stream and its model, respectively. The resulting fitness function / combines the 
above defined components in the following way: 



/ 



(A5) 



Our definition of the fitness function is different from the approach by I Ruzicka et al. ( 20091) . Moreover, 
the additional fitness component has been introduced by Equation (|A3I) . Therefore, the value of the 
fitness function returned for the given model needs further analysis to clarify whether the identified local 
peaks of the fitness function correspond to satisfactory models. This means that every system has its own 
threshold value of the fitness function which needs to be defined. 
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Table 1 

Free parameters of the LMC-SMC-Milky Way interaction 



Parameter 


Value 


Comment 


(CKlmc, <5lmc) 


(81.90° ±0.98°, -69.87° ±0.41°) 


Equatorial coordinates 


(f^smc > <5smc) 


(13.2° ±0.3°, -72.5° ±0.3°) 




(m - M) Imc 


18.5 ±0.1 


Distance moduli 


(m - M) 

\ ' smc 


18.85 ±0.10 




^[masyr- 1 ] 


(-2.11,-1.92) 


Proper motion components 


M^masyr- 1 ] 


(±0.39, ±0.49) 




^[masyr- 1 ] 


(-1.34, -0.69) 




MN^fmasyr- 1 ] 


(-1.35,-0.99) 






262.2 ±3.4 


Line-of-sight velocities 


^Tdkms- 1 


146.0 ± 0.60 




r&pc] 


(9.0,12.0) 


Initial radii of the particle disks 


rSSfkpc] 


(5.0,8.0) 




(^lmc 5 2^1mc) 


(34.7° ±6.2°, 129.9° ±6.0°) 


Inclination and position angles of the particle disks 


(ismc ; Psmc) 


(60° ±20°, p = 45° ±20°) 




Ootkms -1 ] 


(210,260) 


LSR circular velocity 


9 


(0.71,1.30) 


Flattening of the DM halo of the Milky Way 
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Fig. 1. — The 2D projections of the Magellanic parameter space to the (/x N , /z w )-plane for both the LMC 
(left plot) and the SMC. The gray fillings mark the proper motion rang es explored by ge n etic a lgorithm. 

3bJ)(K06a), 



(20061 



The labels indicate the proper motio ns as expect e d by the studies by iKallivavalil et "al 
iKallivavalil et al" (|2006al )fK06b'). and IPiatek et all (|2008)(PI08). The ellipses show the 68.3% confidence 



regions. 
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Fig. 2. — View of the Magellanic System. The po- 
sitions of the Sun and of the Magellanic Clouds are 
shown in a galactocentric Cartesian frame. The 
frame was selected so that its z-axis coincides with 
the axis of the Milky Way disk and points towards 
the Northern Galactic Pole. The current Solar po- 
sition vector is R Q = (— 9 /^o, 0, 0) kpc, and the 
Sun is moving in the direction of the y-axis, i.e. 
v® = (0,6 ,0)kms _1 . The Magellanic Clouds 
are ps20kpc far from each other, and their present 
day velocity vectors are almost parallel. 
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Fig. 3. — Galactocentric velocities of the Magel- 
lanic Clouds. The magnitude of the velocity vector 
of the LMC as a function of the LSR circular veloc- 
ity is displayed in the upper plot. The lower plot 
shows the same relation for the SMC. Both plots 
assume the galactocentric coordinate frame de- 
picted in Figure ((2|) and the L SR angular rotation 
rate lip = 29.45 km s -1 kpc -1 (jReid fc Brunthaler 



2004). Both Clouds slow down with respect to 



the Galactic Center if the LSR circular velocity is 
increased. 
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Fig. 4. — Distribution of all high-fitness models of the Magellanic System identified by the genetic algorithm. 
The upper row shows the fitness of the models as a function of the western (left plot) and of the northern 
LMC proper motion component, respectively. The lower row offers the same relations for the SMC. 
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Fig. 5. — Total mass of the Milky Way enclosed 
within the radius of 250 kpc as a function of the 
LSR circular velocity 0o and of the DM halo flat- 
tening q. All the high-fitness models identified by 
the genetic algorithm search are over-plotted ac- 
cording to their [q, &q] positions. While no trend 
exists regarding the dependence on O , the models 
for oblate (q < 1) haloes are preferred. 
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Fig. 6. — The galactocentric velocity of the LMC 
as a function of the total mass of the Milky Way. 
The models of the interaction between the Mag- 
ellanic Clouds and the Galaxy which were identi- 
fied by the genetic algorithm are divided into four 
groups regarding the value of their fitness. The 
curves (solid lines) corresponding to the paramet- 
ric equations (|T51) and (|T9")) are plotted for several 
randomly selected models of the fitness / > 0.60 
(squares with crosses inside), i.e. for different 
combinations of the parameters q, (c%q C , Si mc ), 
(m - M) lmc , (^ c ,/i^ c ), and v\Z- The dotted 
lines demonstrate the proximity of the parametric 
curves to a linear relation. 
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Fig. 7. — Orbits of the Magellanic Clouds as functions of the Milky Way halo flattening q (left column) and of 
the LSR circular velocity 9o, respectively. The parameter q was varied within the range of (0.71, 1.30) while 
Oo was fixed to the value of 238 km s -1 . When the LSR circular velocity was the variable (8o = (210, 260)), 
the flattening q — 0.80 was selected (right column). The upper row depicts the past time dependence of the 
LMC-Galaxy and the SMC-Galaxy distances. The plots in the lower row show the LMC (more to the right) 
and SMC orbits over the last 300 Myr projected to the plane of sky. Galactic coordinates are used. The 
color coding is such that the black to light gray transition corresponds to changing the actual variable from 
its minimum to its maximum. The c ontour maps s hows the low resolution HI observations of the Magellanic 



Clouds and associated structures bv lBriins et al.l (2005) 
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Fig. 8. — LSR radial velocity profiles along the 
orbital tracks of the Magellanic Clouds. The ra- 
dial velocity is plotted a s a function of the Magel- 
lanic longitude (see, e.g. lBriins et alJl2005j ) for the 
models identified by the genetic algorithm. The 
best models of / > 0.60 (upper plot) and the low- 
quality fits of / < 0.50 are shown, respectively. 
The past positions of the Magellanic Clouds are 
marked at the times of —150 Myr and —450 Myr. 
The bold line stretched between the Magellanic 
longitudes of 0° and —110° corresponds to the ob- 
served LSR radial velocity profile of the Magellanic 
Stream. 
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Fig. 9. — Intensity of the gravitational field of the 
logarithmic halo. The radial (solid line) and the z 
components of the gravitational acceleration due 
to the axially symmetric logarithmic halo of the 
Milky Way are plotted as functions of the halo 
flattening parameter q. The values and the ratio 
of the components of the gravitational acceleration 
depend on the position as well. The radial force 
increases if the halo flattening is increased while 
the axial component of the gravitational field de- 
creases in such a case. 
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Fig. 10. — Fitness of the models identified by our 
evolutionary optimizer is plotted as a function of 
the scalar product of the current LMC and SMC 
velocity vectors. In most cases the deviation is 
smaller than 15°. 
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Fig. 11. — Time dependence of the distance 
between the centers of mass of the Magellanic 
Clouds. The function is plotted for the models of 
/ > 0.55 back to the time of —0.7 Gyr. In all cases 
the LMC-SMC distance is increasing at present 
implying a minimum of the distance occurring in 
the past. 
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Fig. 12. — Current rate at which the distance between the Magellanic Clouds changes with time. The 
function dD\_ s /dt is plotted as a function of the SMC proper motion components /xw and /^n- The proper 
motion components of the LMC correspond to those of a selected high-fidelity (/ > 0.60) model (left plot), 
and to the limiting values of the ranges given by Equation (TTJ), respectively. The minimum value for the 
western proper motion and the maximum of the northern LMC proper motion component were chosen to 
achieve the minimum spatial velocity. At present the rate dD\- s /dt is positive. The distance between the 
Clouds is increasing. 
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Fig. 13. — Time dependence of the distance 
between the centers of mass of the Magellanic 
Clouds. The function is plotted for the models 
of / > 0.60 back to the time of —4.0 Gyr. A close 
(d ~ 15kpc) encounter between the LMC and the 
SMC occurred at t < —2.5 Gyr in all cases de- 
spite the fact that the Clouds cannot be consid- 
ered gravitationally bound to each other until the 
recent past. 
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Fig. 14. — The threshold of the fitness function. The plots on the left hand side present an example of a 
satisfac tory reproduction of the H I observations of the large-scale structures associated with the Magellanic 
Clouds (|Bruns et alj|2005l) . The plots in the right hand column illustrate the typical output from a model for 
the LMC-SMC-Galaxy interaction that is not considered successful. The upper row of this figure offers the 
simulated low-resolution 3D data-cube for both models. The column density isosurface of S = 10 _4 £ max 
is depicted, together with the contour plots of the integrated relative column densities of HI projected 
to the position-position and position-LSR radial velocity spaces, respectively. The lower row provides a 
different view of the position versus LSR radial velocity dependence for the modeled Magellanic Stream. 
Positions of all individual LMC/SMC particles are plotted. Its LSR radial velocity profile is compared to 
the mean profile of the obser ved Stream (thick black line). For the definition of the Magellanic longitude, 
sec IWannier fc Wrixonl fll972h . Positions of the centers of mass of the Clouds are indicated at the time of 
— 150 Myr and —300 Myr, respectively. 
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Fig. 15. — Large-scale distribution of H I associated with the Magellanic Clouds - observed versus simulated. 
The comparison of the observed (left column) low-resolution HI data-cube for the extended Magellanic 
structures (the Magellanic Stream, the Leading Arm) and of its modeled counterpart. The plots in the 
upper row depict the HI column density isosurfaces of £ = 10~ 4 £ max , and the contour plots of the integrated 
relative column densities of HI projected to the position-position and position-LSR radial velocity spaces, 
respectively. The middle row shows the mean column density of HI in the Stream as a function of the 
Magellanic longitude. The third row offers a detailed view of the integrated HI column densities projected 
to the plane of sky. The orbital tracks of the Clouds for the actual model are over-plotted over the last 
300 Myr (lower left plot) and over the last 3Gyr, respectively. The solid line corresponds to the past orbit 
of the LMC, while the dash-dotted line was used for the SMC. The present western and norther proper 
motions for both the Clouds are indicated. 3g 
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Fig. 16.— Relative number of the LMC/SMC 
particles strongly disturbed due to the encounter 
events in the LMC-SMC-Galaxy interaction for 
a high-fitness model of / = 0.61 (see also Fig- 
ure The upper plot shows the relative counts 
of the disturbed particles in eight bins of 500 Myr, 
starting at the time of —4 Gyr. The last bin is 
subdivided into five equal sections. The lower plot 
depicts the present distribution of the LMC/SMC 
particles. They are color-coded according to the 
epoch when they were stripped from the Clouds. 
The particles that remained bound to their galaxy 
of origin over the entire period of the simulation 
were not plotted. The Magellanic Stream is a com- 
position of filaments of different ages spread along 
the past orbits of the Magellanic Clouds. 
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Fig. 17. — Tidal acceleration for an axially sym- 
metric logarithmic potential as a function of its 
flattening parameter q. The values are calculated 
for two points on circular polar orbits of a con- 
stant radial distance Ar = lOkpc. Three orbits of 
the radii of 50kpc, 100 kpc and 150 kpc are shown, 
respectively. The plots show the mean value, the 
maximum and also the minimum of the accelera- 
tion on the given orbits. 
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